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Abstract. The particle acceleration process at a shock 
wave, in the presence of the second-order Fermi accelera- 
tion in the turbulent medium near the shock, is discussed 
as an alternative explanation for the observed flat syn- 
chrotron spectra of supernova remnants (SNRs) in molec- 
ular clouds. We argue that medium Alfven Mach number 
shocks considered by Chevalier (1999) for such SNRs can 
naturally lead to the observed spectral indices. 
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1. Introduction 

In a recent paper Chevalier (1999) discusses physical con- 
ditions in supernova remnants (SNRs) evolving inside 
molecular clouds. He considers a clumpy gas distribution 
in the clouds, where the dominant part of the cloud mass 
is contained in the compact dense clumps filling only 10% 
or less of the cloud volume and the rest of the cloud is 
filled with tenuous gas with a number density N ~ 10 
cm -3 . Thus a SNR exploding inside a cloud evolves mostly 
in such low density medium. Chevalier discusses a num- 
ber of consequences of such a model and compares it to 
observations of three SNRs: W44, IC 433 and 3C391. In 
his discussion the observed, very flat, cosmic ray electron 
distributions responsible for the radio synchrotron spectra 
are interpreted as shock-compressed ambient distributions 
radiating downstream of the radiative shock. 

However, the existence of a very flat ambient electron 
distribution in a molecular cloud is a matter for debate. 
The uniform magnetic field structures observed in several 
clouds and efficient damping of short Alfven waves in a 
partly neutral medium (Hartquist & Morfill 1984) suggests 
a possibility of an efficient cosmic ray exchange between 
clouds and a steep-spectrum galactic cosmic ray popula- 
tion. In the present note we point out the possibility of an 



acceleration of the flat spectrum electrons at a SNR shock 
wave, if the second-order Fermi acceleration in the vicin- 
ity of the shock is taken into account. In the next section 
we summarize results of an approximate analytic theory 
for the particle spectral indices (Ostrowski & Schlickeiser 
1993, = 'OS93'), with a few typing errors of the original 
paper corrected. Then, in section 3, we demonstrate that 
the physical conditions considered by Chevalier (199) for 
shock waves - with the Alfven velocity, Va, non-negligible 
in comparison to the shock velocity - allow for genera- 
tion of the required flat particle distributions (c.f, also, 
Drury 1983, Hartquist & Morfill 1983, Droge et al. 1987; 
Schlickeiser & Fifrst 1989). In the final section (Section 4) 
we briefly discuss an application of the present model to 
actual conditions in astrophysical objects. 



2. Derivation of the particle spectral index 

OS93 derived a simplified kinetic equation for the parti- 
cle distribution formed at the parallel shock wave due to 
action of the first-order acceleration at the shock and the 
second-order acceleration in the shock turbulent vicinity. 
With the phase-space distribution function at the shock, 
f(p) = f(p,x = 0), an 'integral' distribution is defined as 



F(p) = f(p) 



Klip) , «2(P) 



Ul 



(2.1) 



where Ui and Ki are the plasma flow velocity in the shock 
frame and the spatial diffusion coefficient of cosmic ray 
particles, respectively (i = 1 upstream and 2 downstream 
of the shock) . Analogously, the momentum diffusion coef- 
ficient is indicated as Di (i = 1, 2). For the simple case of 
Ki — const (i — 1,2), and D\ = Di = D the approximate 
transport equation for the function (2.1) takes the closed 
form: 
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p 2 dp 



p 2 D(p)-F(p) 



l__d 

p 2 dp 



r(p) 



(2.2) 



where we have included the source term Q(p). The mean 
acceleration speed due to the first order process at the 
shock (-^f) and the mean escape time due to advection 
downstream of the shock r(p) can be derived from the 
spatial diffusion equation (cf. OS93) as: 



R-l 
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where the shock compression R = U1/U2. In more gen- 
eral conditions with (p) the function f(p) must be 
explicitly present in the kinetic equation: 
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Let us note in the above equation the asymmetry between 
the first two terms. If only the first-order acceleration 
takes place ( Di = = D 2 ), the known solution /(p) oc p a , 
with a = -3R/iR - I), is reproduced by Eq. (2.5). Ob- 
taining a solution for a more general situation may be 
a difficult task. However, if we consider momenta much 
above (or below) the injection momentum, and the power- 
law form for the diffusion coefficient n oc p v , the solution 
is also a power-law. In such conditions, from Eq. (2.5) one 
can derive the spectral index a. The Skilling (1975) for- 
mula is used for the momentum diffusion coefficient, re- 
lating it to the spatial diffusion, Dip) = V^p 2 / (9k(p)) . 
For a given jump condition at the shock for the Alfvcn ve- 
locity and a given value of 77, the resulting spectral index 
depends on only two parameters, the shock compression 
ratio R and the velocity ratio Va,i/Ui. 

OS93 checked the validity range of the approximate 
equation (2.5) with the use of numerical simulations. For 
77 7^ the equation can be used for the range of parameters 
preserving Va « U, while, for k = const, it provides a 
quite reasonable description of the particle spectrum at 
all V A < U. 

In order to find a power-law solution of Eq. (2.5) let us 
assume the following forms for the diffusion coefficients (i 



= 1, 2): Kiip) = Kojp 11 , Diip) = D .ip 2 11 , where the con- 
stants Ko,i and Z?o.i are related according to the formula 
Do,i — ^4,i/(9 K o,i)- With these formulae and the power- 
law form for the distribution function, fip) = f p a , equa- 
tion (2.5) yields a quadratic equation for a : 



aa + ba + c = 



with coefficients: 



a = — 
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Equation (2.6) has two solutions valid for, respectively, 
particle momentum much below (<r > 0) and much above 
(a < 0) the injection momentum. Only the later one is of 
interest in the present considerations. Let us also note that 
the energy spectral index often appearing in the literature 
is F = a — 2 and the synchrotron spectral index a = 
(a-3)/2=(r-l)/2. 

3. Spectral indices in realistic SNRs 

In the SNR model discussed by Chevalier (1999) the shock 
wave propagates inside an inhomogeneous cloud. In the 
cloud, the dense clumps occupy less than 10% of the 
cloud volume and the shock evolution proceeds mostly in 
the tenuous inter-clump medium with a number density 
N ~ 10 cm" 3 and a magnetic field B w 2 • 10~ 5 G. Much 
stronger magnetic fields (~ 10~ 4 G) may occur further 
downstream in the radiative shock. The measured shock 
velocities are in the range Ui s=s 80 - 150 km/s. With 
the notation N = A/(10cm- 3 ), B = B/{2 ■ 10~ 5 G) and 
Uq = C/i/(100km/s) we find the ratio 



Va^ 
Ui 



= 0.UB N o 1/2 C/ ( 



(3-1) 



It can be of order 0.1 for the parameters considered above. 
Thus, the second order acceleration process can substan- 
tially modify the energy spectrum of particles accelerated 
at the shock. In Fig-s 1,2 we present spectral indices de- 
rived from the formulae of the previous section for the 
'canonical' Kolmogorov value n = 0.67 . The shock com- 
pression ratio R = 4.0 for the high Mach number adia- 
batic shock is assumed. An efficient particle acceleration 
requires a high amplitude turbulence near the shock (cf. 
Ostrowski 1994) and oblique magnetic field configurations 
may occur there. Thus, the mean magnetic field down- 
stream of the shock is B 2 > B x due to shock compression 
of both the uniform and the turbulent field components. 
Below we use an effective jump condition for the Alfvcn 
velocity Va,i = Va,i (in general B\ < B 2 < B\R and, 
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respectively, Va,i/VR < Va,2 < Va,iVR)- The results for 
Va,i/Ui < 0.2 are considered, where the analytic formulae 
provide an accurate approximation for the actual spectra. 
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Fig. 1. A map of the synchrotron spectral index a in 
co-ordinates of the magnetic field B given in fiG and the gas 
number density N in cm~ 3 , both taken upstream of the shock. 

For mean models for the SNRs discussed by Chevalier 
(1999) one has: 

a. ) for W44: N a = 0.5, U = 1.5, a = 0.33 (r = 1.72) 

b. ) for IC 443 (shell A): N = 1-5, U = 0.8, a = 0.36 
(r = 1.66) 

c. ) for 3C391: N = 1.0, U = 3.0, a = 0.55 (T = 2.1) 
With these parameters one can easily explain flat spectral 
indices for the cases (a) and (b) if reasonable values of 
the magnetic field are involved. Fig. 2 shows the spectral 
index a versus the magnetic field B for the above listed 
choices of particle density N — N and the shock veloc- 
ity Ui = Uo- The measured indices are indicated on the 
respective curves. The steep spectrum of 3C391 can not 
be exactly reproduced with the use of the simple model 
considered in this paper, as even neglecting the second or- 
der acceleration leads, for the assumed shock compression 
R = 4.0, to a somewhat flatter spectrum with a — 0.5. 
However, the observed trend in the spectral index varia- 
tion is well reproduced by the model. 

4. Discussion 

A derivation of the particle spectral index at a shock front 
is presented for a situation involving the second-order ac- 
celeration process in the shock vicinity. We prove that 
even very weak shocks may produce very flat cosmic ray 
particle spectra in the presence of momentum diffusion. 
As a consequence, the dependence of the particle spectral 
index of the shock compression ratio can be weaker than 
that predicted for the case of pure first-order acceleration. 
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Fig. 2. The synchrotron spectral index a versus the magnetic 
field B for the respective values of N and Ui. The curves for 
different objects are indicated with the respective SNR symbol. 
The observed spectral indices are indicated on the curves for 
IC 443 and W44. 

One should note an important feature of the acceler- 
ation process in the presence of the second-order Fermi 
process acting near the shock. Because the same Alfven 
waves scatter particle momentum in direction and in mag- 
nitude, there exists a strict link between the first- and 
the second-order acceleration processes. The particle spec- 
trum is shaped depending on Va,i/U\, by the compression 
R, the momentum dependence of k and, possibly, by the 
anisotropy of Alfven waves determining the ratio of k to 
D. The last factor is not discussed here since under the 
present considerations we restrict ourselves to isotropic 
wave fields, but the role of anisotropy consists in decreas- 
ing the importance of momentum diffusion. Only through 
the above parameters the spectrum can be influenced by 
other physical characteristics of the shock and the medium 
in which it propagates. It is important that the Alfven 
waves' amplitude, or the magnitude of the spatial diffu- 
sion coefficient related to it, do not play any substantial 
role in determining the spectral inclination as both, first- 
and second-order, processes scale with it in the same way. 
One could argue that an energy transfer from waves to 
particles can cause quick damping of the waves and will 
leave us with the pure first-order Fermi acceleration at 
the shock. From the above discussion one can infer that 
such an objection is not valid for the isotropic wave field. 
Only the presence of high amplitude one-directional wave 
field enables the efficient first-order shock acceleration in 
the absence of momentum diffusion, but a detailed cal- 
culation is not possible without a detailed model for the 
upstream and downstream wave fields. 

In comparison to parallel shocks, the magnetic field in- 
clined to the shock normal leads to a higher mean energy 
gain of particles interacting with the shock and a higher 
escape probability for downstream particles. Also, due to 
small cross-field diffusion the normal diffusive length scale 
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near the shock decreases. As long as the shock veloc- 
ity along the field is non-relativistic, the particle spec- 
trum produced in the first-order process is not influenced 
by a field inclination. Because of smaller diffusive zones 
near the shock, the role of momentum diffusion may be 
of lesser importance as compared to the parallel case. 
This effect can be partly weakened by the presence of the 
magnetic field compression at the shock, which leads to 
higher Alfvcn velocity downstream of the shock. Also the 
presence of high amplitude Alfven waves (cf. Michalek & 
Ostrowski 1996) and an admixture of fast mode waves 
propagating obliquely with respect to the magnetic field 
(Schlickeiser & Miller 1998, Michalek et al. 1998) may lead 
to a more efficient acceleration than the one considered by 
OS93. Finaly, we would like to note a novel approach to 
the shock acceleration by Vainio & Schlickeiser (1999), 
who discuss conditions allowing for the shock generated 
flat spectra without action of the second order process. 
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